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Simple Lie algebras and topological ODEs 

Marco Bertola, Boris Dubrovin, Di Yang 


Abstract 

For a simple Lie algebra g we define a system of linear ODEs with polynomial coefficients, 
which we call the topological equation of g-type. The dimension of the space of solutions 
regular at infinity is equal to the rank of the Lie algebra. For the simplest example g = SI 2 (C) 
the regular solution can be expressed via products of Airy functions and their derivatives; 
this matrix valued function was used in our previous work [4] for computing logarithmic 
derivatives of the Witten-Kontsevich tau-function. For an arbitrary simple Lie algebra 
we construct a basis in the space of regular solutions to the topological equation called 
generalized Airy resolvents. We also outline applications of the generalized Airy resolvents 
to computing the Witten and Fan-Jarvis-Ruan invariants of the Deligne-Mumford moduli 
spaces of stable algebraic curves. 

Keywords, simple Lie algebra; topological equation; dual topological equation; generalized Airy resol¬ 
vents; normal form; essential series. 


1 Introduction 

In this paper to an arbitrary simple Lie algebra we will associate two families of special functions 
defined by certain systems of linear ODEs with polynomial coefficients along with appropriate 
asymptotic conditions. Let us begin with the first family. 

1.1 Topological ODE and generalized Airy resolvents 

Let (0, [•,•]) be a simple Lie algebra over C of rank n. We sometimes refer to elements of g as 
matrices, even if the constructions will not depend on the choice of a matrix realization of the 
Lie algebra. Choose i) a Cartan subalgebra of g, 11 = {ai,. .. ,an} (Z i)* a. set of simple roots, 
and Ac!)* the root system. Then g admits the root space decomposition 

g = f)©0ga. (1.1) 

«ga 

Let 6 be the highest root w.r.t. 11. Denote by (-I-) : g x g —)■ C the normalized Cartan-Killing 
form [24] such that {0\9) = 2. For any root a, let Ha G 1} denote the root vector of a w.r.t. (-I-)- 
Choose a set of Weyl generators Ei G go,., Fi G g-a^. Hi = 2Hai/{ai\ai) satisfying 

[Ei, Fi] = Hi, [Hi, Ej] = Aij Ej, [Hi, Fj] = - Aij Fj (1.2) 
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where (Aj) denotes the Cartan matrix associated to g. Here and below, free Latin indices take 
integer values from 1 to n unless otherwise indicated. Choose E_e G g-e, Eq G g^- They can 
be normalized by the conditions {Eq \ E_q) = 1 and uj{E_g) = —Eg, where cu : g —> g is the 
Chevalley involution. Let /+ ;= a principal nilpotent element of g. Dehne 

A = /_|_ + A E—g. (T3) 

Consider the following differential equation for a g-valued function M = M(A) of an independent 
variable A 

M'=|Af,A]. ' = T, (1.4) 

Definition 1.1. We call eq. (1-4) the topological differential equation (in short: topological ODE) 
ofg-type. 


Solutions of the differential equation (1.4) are entire functions of the complex variable A. 
At A = cxo they may have a singularity. A solution M(A) is called regular if it grows at most 
polynomially at |A| —)■ oo within a certain chosen sector of the complex A-plane. Denote by 
® 00 ^( 0 ) fhe vector space of regular solutions to (1.4). Our goal is to describe the space ©^®(g) 
for an arbitrary simple Lie algebra g. Actually, we will only deal with asymptotic expansions 
of regular solutions without entering into details about the region of their validity that can be 
described by means of standard techniques of asymptotic analysis of ODEs, see e.g. [35]. 


As a way of introductory example we consider now the simplest instance g = sUiE.). In this 

0 M Write M(A) - ( 

A 0 r ^ ’ I c(A) -a(A) 


case A = 


The equation (1.4) reads 



i / 

V 

= 

\ f 1 

1 \ 

\d 

\ 



(1.5) 


The type of behaviour at inhnity is essentially specified by the eigenvalues of the matrix of 
coefficients. In this case they are equal to 0 and ±2\/A. Solutions corresponding to the nonzero 
eigenvalues have an essential singularity ~ at inhnity. The regular solution is unique 

up to a constant factor. It corresponds to the eigenvector with eigenvalue zero. The solution has 
the form 

A"5 

M(A) = ^ 


This matrix-valued function appeared in our paper [4]. It was used, in a slightly modihed 
normalization, as key tool within an efficient algorithm for computation of intersection numbers 
of tautological classes on the Deligne-Mumford moduli spaces, see eq. (5.14) below. 

The topological ODE for the Lie algebra s/ 2 (C) is closely related to the theory of Airy func¬ 
tions. Indeed, from (1.5) it readily follows that 

a = ^6', c = -^b" + Xb (1.7) 


/ (6g-5)!! ■> 3q+2 r) y^oo (Gg-p!! \-3g \ 

Z^g=l 969-1.(g-i)!^ A^g=0 969-g! ^ 

y^oo 6g+l (6g-l)!! \-3q+l _ y^°° (6g-5)!! \-3q+2 

\ ^ Z^g=0 6g-l 969-g! ^ 2.^g=l 969-1-(g-l)!^ y 


( 1 . 6 ) 
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while for b = b{X) one arrives at a third order ODE 


- h'" + 2Xb' + b = 0. 

Solutions to this equation are products of solutions to the Airy equation 


( 1 , 8 ) 




(1.9) 


That is, b{X) is the diagonal value of the resolvent of the Airy equation. The matrix M(A) can 
also be considered as resolvent of the matrix version of the Airy equation 


■ d / 

' 0 

1 M 

.dA + ( 

.A 

o 

1_ 



( 1 . 10 ) 


Let us proceed now with the general case. The eigenvectors of the matrix of coefficients of 
the system (1.4) with zero eigenvalue correspond to matrices commuting with A. As observed in 
[27], the kernel of the linear operator 


adA(Ao) : 0 0 

for any Aq G C has dimension n = rkg. This suggests that the dimension of the space of regular 
solutions to the topological ODE (1.4) is equal to n. We show in Thm. 1.2 that this is the case 
indeed. Moreover, we will construct a particular basis in the space ©^®(g). To this end denote 
-^( 0 ) = 0 ® C[A, A“^] the loop algebra of 0 . Introduce the principal gradation on the Lie algebra 
L(g) in the following way 

degEj = l, degi7j = 0, degFj =-1, i = 1,..., n, degA = h. (1-11) 

Here h is the Coxeter number of g. According to this gradation we have deg A = 1. Recall [23] 
that KeradA C L(g) has the following decomposition 

n 

Ker adA = © CAj, deg Aj = j & E ■= |_|(mj + hlX) (1-12) 

j&E i=l 


where the integers 


1 = mi < 1712 < • • • < 'mn-l < = h — 1 (1-13) 

are the exponents of g. The matrices A* commute pairwise 

[Ai,Aj] = 0, yiJeE. (1.14) 

They can be normalized in such a way [25, 24] that 

X^ma+kh = JXma ^ fc G Z, (1-15) 


In particular, we choose Ai = A. 
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Theorem 1.2. 1) dime = rkg. 

2) There exists a unique basis Mi{X), .. ., Mn{X) in the space of the form 


Ma{X) = \ 


mg 

h 


K 


k=l 


A'4,i(A) 6 i(g), deg M„,t(A) = m„ - (A + 1)*:, a=l,...,ij. 


(1.17) 


The proof is given in Sect. 2. 

Definition 1.3. The solutions (1.17) to the topological ODE of Q-type are called generalized Airy 
resolvents of Q-type. 


The topological ODE (1.4) can be recast into the form 


A. 

dX 


+ A,M 


0 . 


(1.18) 


Using such a representation along with (1.14), (1.16) as well as the invariance property of the 
Cartan-Killing form one arrives at 

Proposition 1.4. The generalized Airy resolvents commute pairwise 


[Ma{X),Mb{X)] = 0, a,b=l,...,n. 


(1.19) 


Moreover, {Ma{X) \ Mb{X)) = h 6 a+b,n+i- 

Example 1.5. Consider the A^ case, g = s/„+i(C). The (n + 1) x (n + 1) matrix A reads 


A 


/ 0 1 0 

0 0 1 


0 0 ■■ 

\ A 0 0 


0 \ 
0 


1 

0 / 


The exponents are ma = a, a = 1,... ,n, the Coxeter number is h = n 1. One can choose 


Then 


Aj = A*, i = 1,... ,n. 


AIa{X) = X "+l 


A“ + 


(n + 1) A ^ 4 

^ ' V=1 


^ 2-va _ 1 


sin 


2 TTU ^ 

n+l 


( 1 . 20 ) 


where ( = e^ is a root of unity, = diag (1, ..., C”^) • The recursion procedure for 

computing subsequent terms will be explained below, see Example 2.4- In the particular case n = 1 
one obtains the solution (1.6). 
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The recursive procedure for computing the series expansions of the generalized Airy resolvents 
outlined in the Example 1.5 above can be represented in an alternative form. The basic idea 
can be already seen from the sl 2 {C) example: the topological ODE for a matrix-valued function 
reduces to a scalar differential equation (1.8) for one of the entries of the matrix (called b{X) 
above); other matrix entries are expressed in terms of b{X) and its derivatives, see eq. (1.7). 

Let us explain a general method for reducing the topological ODE for an arbitrary simple Lie 
algebra of rank u to a system of ODEs with n dependent variables. To this end we will use the 
Jordan decomposition of the operator ad/^ described hereafter. 

Let G 1) be the Weyl co-vector of g, whose dehning equations are 

ai{p^) = 1, i = 1,... ,n. (1-21) 

Write = Ym=i Hi, Xi G C and dehne J_ = 2 Ym=i Hi. Then /+,/_, form an s/ 2 (C) Lie 
algebra: 

[p\ /+] = /+, [p\ /_] = [/+, /_] = 2p\ ( 1 . 22 ) 

According to [27, 2], there exist elements 7 ^,..., 7 " G g such that 

Kerad/_ = Spanc{ 7 \ ... , 7 "}, [p'^, 7 *] =- 771 ^ 7 *. (1.23) 

Fix { 7 ^,..., 7 "}, then the lowest weight decomposition of g has the form 

n 

B = 0r. r = Spant.{T‘.ad,,T’.....ad77}. (1.24) 

i=l 


Here each D is an s/ 2 (C)-module. Any g-valued function M(A) can be uniquely represented in 
the form 

n n 2mi—l 

M(A) = + E E ad" 7* (1.25) 

2=1 2=1 m =0 

where 5'j(A), Ki^i^X) are certain complex-valued functions. 

Theorem 1.6. Let M he any solution of (1.4)■ The functions Kim have the following expressions 


n 2mi — l—m 


n mi — l—m 


u=l v=0 u=l v=0 

^ (-1)" m ^ m,...., 2m, - 1 

(1.26) 


with constant coefficients k^muvi ^imuv independent of the choice of the solution M. Moreover, the 
topological ODE (1-4) is eguivalent to a system of linear ODEs for Si,..., Sn of the form 


^(2?72i + l) 


n 2mi 

EE 

u=l 2;=0 


kL sp + feL 


sp, i = l, 


U=1 V=1 


Here are constants independent of the choice of the solution M. 


(1.27) 
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The proof is given in Sect. 3. Note that the coefficients coincide with k^rnuvi Hmuv 

m = — 1 . See details in the proof. 

We call eqs. (1.27) the reduced topological ODEs of g-type. Denote by S'a;i(A),..., S'a;„(A) the 
S'-coefficients in the decomposition (1.25) of the generalized Airy resolvent Ma{\), a = 1, ... ,n. 

Definition 1.7. The series Sa-,i{X) are called essential series of Q-type. 

Example 1.8 (The Ai case). The essential series Ai;! coincides with the {l,2)-entry of the Airy 
resolvent (1.6) (the function b{X) in (1.5)). The expression (1.26) is given explicitly by (1.7). 


Example 1.9 (The A 2 case). The decomposition (1.25) reads 

/ 27711 + K22 -2^1 - 37723 6^2 \ 

M = 277io + 7721 -27722 37723 - 2^i . (1.28) 

\ K 20 277io — K 21 — 277ii + K 22 ) 

The expression (1.26) is given explicitly by 

77io = 5" + 3A^2, Kn = -S[, (1.29) 

772o = ^f-2AAi, K2i = -S”', K22 = S'), K23 = -S',. (1.30) 

The reduced eguations (1.27) read 

6^2 + 2^f^+ 9A^2 = 0, (1.31) 

2^1 + 2^J^^ + 15A' - A® + A (9^'' + 6S[) = 0. (1.32) 


Explicit expressions of the essential series Sa-j of the A 2 -type can be found in Sect. 5. 
Proposition 1.10. 1. Eor any a G {1,2,.. .n], the essential series Sa-i satisfy 

1) ^,;i(A),...,^,;„(A)eA-"^-C[[A-i]]; 

2 ) the vectors {Sa.,i)i=i...n, a = 1,... ,n are linearly independent] 

3) Sa-a is non-zero. 

2. Let 0 Dn= 2 k and M be a regular solution of the topological ODE (I. 4 ) such that the Sa- 
coefficient of M coincides with Sa-a,' then M = Ma. If g is of type Dn with n even, the statement 
remains valid for a ^ {n/2, n /2 + 1 }, and it should be modified for a = n/2 and a = n /2 + 1 in 
the following way: M = Ma if the Sn/ 2 , Sn/ 2 + 1 -coefficients of M coincide with Sa-n/ 2 , Sa.,n/ 2 +i, 
respectively. 

The proof will be given in Sect. 3. We call a = 1,..., n the fundamental series of g-type. 

1.2 Dual topological ODE and its normal form 

Solutions to (1.4) can be expressed by a suitable integral transform via another class of func¬ 
tions satisfying somewhat simpler differential equations. These functions will be constructed as 
solutions to another system of linear ODEs for a g-valued function G = G{x), a; G C. In the 
notations of the previous section it reads 

\G', E—q\ -\- \G, 7+] -\- X G = D, 
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(1.33) 


Definition 1.11. We call equation (1.33) the dual topological differential equation (dual topo¬ 
logical ODE for short). The space of solutions of (1.33) is denoted by G{q). 

Solutions to the dual topological ODE are related with the generalized Airy resolvents by the 
following Laplace-Borel transform 

M(A) = [ G{x)e-^^dx, (1.34) 

Jc 

where O is a suitable contour on the complex x-plane, which can depend on the choice of a 
solution M(A) G 

Example 1.12. For g = sOiC) the dual topological equation for the matrix-valued function 

G = ( “hi ) reduces to 

\ 3W) -a{x) ) 


a'= 



(1.35) 


This yields the following integral representation^ for solutions to the differential equation (1.8) 

OO 

6 (A) = ^ (1.36) 

Vvr J \Jx 

0 


For an arbitrary simple Lie algebra we will now describe the reduction procedure for the dual 
topological ODE similar to the one used in the previous subsection. Write 

n n 2mi—1 

G{x) = Y, 0,(x) ad^^y + (1-37) 

2=1 2=1 m =0 


(cf. (1.25) above). 

Theorem 1.13. i). Let G be any solution of (1.33). Then Kim have the following expressions 

n 2mi—m n rrii—l—m 

Kim = Y ^irnuv Y ^ 2:''“ V«) , m = 0,..., m* - 1, 

u=l V=0 U=1 21=0 

Kim = y, m = ruj,..., 2mi (1.38) 


^ Here and below the Laplace-type integrals are understood in a formal sense by applying a term-by-term 
integration 

1 


OO 


b1_xx dx _ 1 


_ eTT--- _= _ V / - f —^ e-^^— = — V y 3fc_j 

sJL sJT yj k\\vi) sJL klU>^ 

Q K -U Q K -U 


In order to arrive at an integral representation of a solution with the same asymptotic expansion within a certain 

{ 2 fc 4 -l) 7 ri 

sector near A = oo one can integrate over a ray of the form x = r e 3 ,0<r<oo for /c = 0, 1, or 2. 
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where k^rnuv) Hmuv same constants as in Thm. 1.6. Moreover, the dual topological ODE 

(1.33) is equivalent to a system 0 /first order ODEs 


n 2mi 


mi 


X 






•V kO 


u=l v=0 


^iuv ^ 


^inv a:’' Vn + (j)'J = 0. 


u=l V = 1 


(1.39) 


Here, i = 1,... ,n, and are constants (the same as in Thm. 1.6). 

a). The point x = Q is a Euchsian singularity of the ODE system (1.39). More precisely, 
equations (1.39) have the following normal form 




X 


2h-2 


+ ^X^Vk(j), (f = {(fi,. . . 


k=0 


1/_1 = diag{-mn+i-a/h)a=i,...,n, Vk e Mat(n, C), k>0. 


(1.40) 

(1.41) 


Hi). dimc^(0) = rkg. 

Remark 1.14. Choosing the contour C in several ways in (1.34) one obtains a complete basis 
of solutions of ( 1 . 4 )- 

Example 1.15. Eor g = sOiC) one reads from (1.35) that V-i = — Vq = Vi = 0, V 2 = j. 
Example 1.16. Eor g = s/ 3 (C) the normal form (I. 40 ) of the dual topological equation reads 


9V20J+6I 



(1.42) 


Let Ga denote the Laplace-Borel transform of Ma, and 0a;i,..., 4>a-,n the 0-coefficients of Ga, 
a = 1,... ,n. Then the matrix $ defined by <hja := 4>a-,i is a fundamental solution matrix of the 
normal form (1.40). 

Definition 1.17. We call <pa-i the dual essential series, and c()a-,a the dual fundamental series. 

Conjecture 1.18. Eor any a E {l,...,n}, the dual fundamental series (fa-,a satisfies a scalar 
linear ODE of order less than or equal to n with polynomial coefficients in x. Such an ODE with 
the smallest possible order, say ra, is called a dominant equation and ra is called a dominant 
number. The point x = 0 is a regular singularity of the dominant equation. If h is a prime 
number, then ri = ... = = h — 1. 

For low ranks, the above conjecture is verified, and the dual fundamental series (fa-,a can be 
expressed via elementary functions or combinations of elementary and Bessel functions, see Table 
1. However, already for the case of the Lie algebra s/ 5 (C) we were not able to identify the dual 
fundamental series with a classical special function. 

Last but not least: why do we call eg. (I. 4 ) the topological differential equation? For the 
case of g = sOiC) the Airy resolvent (1.6) appears in the expression of [4] for generating series 
of intersection numbers of "^-classes on the Deligne-Mumford moduli spaces A4g j\f; see also the 
formula (5.14) below. It turns out that for a simply-laced Lie algebra g a representation similar 





Type of g 

Special functions arising in dual fundamental series 

Ai 

exponential function 

rank 2 

Bessel functions 

rank 3 

Bessel functions, exponential functions 

A 4 

solutions to scalar ODEs of order 4 

D 4 

Bessel functions, exponential functions 


Table 1: Special functions arising from the dual topological ODE of g-type. 

to that of [4] can be used for computing the E. Witten and H. Fan-T. Jarvis-Y. Ruan (FJRW) 
intersection numbers. We give the precise expressions in the last section of the present paper, 
see Theorems 5.1 and 5.5 below. Connection between FJRW invariants and tan-functions of the 
Drinfeld-Sokolov hierarchy oi AD F-type conjectured in [38] and proven in [16] and [17] plays an 
important role in derivation of these expressions. The non-simply-laced analogue of the FJRW 
intersection numbers, according to the recent paper by S.-Q. Liu, Y. Ruan and Y. Zhang [30], is 
also related to tan-functions of the Drinfeld-Sokolov hierarchies oi B C F G-type. 


Organization of the paper. In Sect. 2 we prove Theorem 1.2. In Sect. 3 we prove Thm. 1.6, 
Prop. 1.10 and Thm. 1.13. In Sect. 4, we calculate the dual essential series, in particular the dual 
fundamental series for several special cases. Concluding remarks are given in Sect. 5. We give 
definition and examples of generalized Airy functions in Appendix A. 
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Hamiltonian systems in hnite and inhnite dimensions” of Italian Ministry of Universities and 
Researches. M. B. is snpported in part by the Natural Sciences and Engineering Research Council 
of Canada (NSERC) and the Fonds de recherche Nature et technologie du Quebec (FQRNT). 


2 Proof of Theorem 1.2 


Before proceeding with the proof, we introduce the gradation operator 

gr = h -F adpv. (2.1) 

Recall that G I) is the Weyl co-vector. The following statement [25] will be useful in the proof. 
Lemma 2.1. An element a & L{q) is homogeneous of principal degree k iff 

gr a = ka. (2.2) 


The action of the gradation operator can be obviously extended to expressions containing 
fractional powers of A. 

The following lemma will also be useful. 
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Lemma 2.2. The following formulce hold true: 


L{q) = KeradA © ImadA, ImadA = (KeradA)'*', (2.3) 

where the orthogonality is with respect to the Cartan-Killing form. 

The proof can be found e.g. in [10]. 

Proof of Thm. 1.2. Let us fix an exponent rua for some a G {1,..., n}. We first prove existence 
and uniqueness of formal solutions of (1.4) of the following form 

OO 

M(A) = A-"^5^Mfc(A), Mo{X) = Am. (2.4) 

k=0 

Mfc(A) G L( 0 ), deg Mk{\) = rUa- {h + l)k (2.5) 

Substituting (2.4) in (1.4) and equating the terms of equal principal degrees we obtain 

[A,Mo] = 0, (2.6) 

7n 1 

[A, Mfc+i] = -M' + -YjMk, k>0. (2.7) 

The hrst condition follows due to the choice Mq = A^^ (see eq. (1.14) above). In order to resolve 
eq. (2.7) for /c = 0 one has to check that 

h A - maAm, G Im adA. 

Indeed, since deg A^^ = a we have 

gr Am, = hX A'^^ + [p^, Am,] = 

So, the equation (2.7) for A; = 0 becomes 

[A, Ml] = ^ [p'^,A^J . 

Due to Lemma 2.2 it suffices to check that the r.h.s. is orthogonal to KeradA. Indeed, 

([p"", I Aj) = {p'^ I [A™,,Aj]) = 0. 

This completes the hrst step of the recursive procedure. 

Assume, by induction, that the k-th term Mfc(A) exists and, moreover, it satishes 

7Tl 1 

-M'f, + -^-Alk G ImadA. 
h X 

Then there exists a matrix M^+i satisfying eq. (2.7). It is determined uniquely modulo an 
element in KeradA. Write 


Mfc+i = A + B, A G ImadA, B G KeradA, 
deg A = deg B = rUa — {h + l){k + 1). 
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( 2 . 8 ) 

(2,9) 


Since the map 


adA : ImadA —)■ ImadA 

is an isomorphism it remains to prove that the matrix B = B{X) is uniquely determined by the 
condition 

Tfl 1 

e ImadA- 

Like above, using the degree condition (2.9) along with the gradation operator (2.1) recast the 
last condition into the form 

[p'^,A + B] + {h + l){k + l){A +B) ± KeradA- (2.10) 

By assumption A _L KeradA- Also [p^, B] ± KeradA due to commutativity of KeradA- Thus eq- 
(2.10) reduces to the system 

{h + l){k + 1) {B I A,) = - ([p^ A] I A,), jeE. (2.11) 

Decompose now B as follows: 

n 

6=1 

Choosing j = mi ,..., and using the normalization (1.16) one hnally obtains 

1 " 

^ ^ ~ Xh{h+l){k + l) ^ I ^-n+l-6) ^rn,. 

This concludes the inductive step and thus completes the proof of existence and uniqueness of 
solutions (1.17) to the topological ODE (1.4). 

The series Mi(A),..., M„(A) are linearly independent since their leading terms are. It is 
known from [35] that there exist analytic solutions to the topological ODE whose asymptotic 
expansions in certain sector coincide with Mj. These solutions we denote again by Mj. We 
conclude that 


Spanc{Mi,...,Mjc6^s(0). (2.12) 

In particular dime ©^^(g) > n. 

On another hand, according to Kostant [27] for any Aq G C the kernel of adA(Ao) : 0 —)■ 0 is of 
dimension n. This implies that 


dime < n. (2-13) 

Hence dime ©]]^^( 0 ) = n and ©]]^^( 0 ) = Spane{Mi,..., Mn). Thm. 1.2 is proved. □ 

The above proof gives an algorithm for computing M^, a = 1,... ,n recursively. 

Remark 2.3. Thm. 1.2 and Prop. I .4 tell that ©[[^^( 0 ) forms a one-parameter family of Abelian 
subalgebras of g. 
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Example 2.4 (The An case). Take the same matrix realization ofg as in Example 1.5, so that the 
normalized Cartan-Killing form coincides with the matrix trace. Recall that h = n + 1, ma = a. 
We have 


Aj = A^, p'^ = diag 


n n — 2 n 




(2.14) 

(2.15) 


2 ’ 2 

i) = Spanc{i/j}"=i. 

Fix an exponent ma. Consider [L, M] = 0, L = dx + A and write M = then 

Mk = Ak + Bk, Ak G ImadA, Bk G KeradA, /c > 0. (2-16) 

Here Ak and Bk have the principal degree degA^ = degi?fc = a — {h + l)k, and 

-1 


Bh — 


^ ^ (adpv^d/j I A^_|_]^_f,) Afj, /c ^ 1. 


A k h{h + 1) 

It can he easily seen that 

Im adA = Spanc[A,A-i] [Hi A^ \ i = 1,... ,n, j = 0,... ,n] , 


Tr(XA^[p\A,]) = 


-h A Xj dn-\-lJ-\-£ j 1) • • • ) ^1 

0 J = 0, 


for X = y^ XjHj. 


(2.17) 

(2.18) 
(2.19) 


2=1 


Thus if Ak (k > 1) is of form Ak = J2j=o^kH^) = I2^=i ^k} (0 < J < n), then 


' '' j=i 

Denote by ( = exp(27r\/—l/h) a root of unity. The map adA acts on ImadA as follows 

adAiXA^) = XX+\ J = 0, ..., n, X,X el), 




/ 2 -1 0 ... 0 \ 


^ Xi ^ 


1 1 -1 ... 0 


^ Xi ^ 

X 2 

— - 

1 0 1 -1 ... 0 


X 2 

\Xn ) 


1 0 ... 1-1 


\Xn ) 


\ 


1 0 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


0 1 / 


Define = diag j-i n+i ’ ^ ■ ■'a, then we have 

adA^n.X) = {C-l)n,X+\ J = d,...,n. 
A direct computation gives 


adpvA° = 0, adpvA^ = [ ^ 


i/=i 


_ f-uil+l) 

C-1 


pu _ ru(e+l) 

adX^adpvA^ = I 1 ^ ^ 1- 


(2.23) 


A^ i> 1; (2.24) 


(2.25) 


(c-ir 

Using these formulce and consider the leading term of M given by Mq = Aa, one immediately 
obtains the formula of the first order approximation (1.20). 
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3 Essential series and dual essential series 


In this section we apply a reduction approach and the technique of Laplace-Borel transform to 
prove Thm. 1.6, Prop. 1.10 and Thm. 1.13. 

First of all, we introduce the lowest weight structure constants. 

Definition 3.1. Fix a choice 0 / 7 ^,..., 7 " and /+ as above. The coefficients uniquely 
determined by 



(3.1) 


are called the lowest weight structure constants of q. 

In this paper we will only use part of these structure constants, namely := The 

above dehnition (3.1) becomes 


n 


2m j 



(3.2) 


Lemma 3.2. The constants are zero unless 

tjb 

s — mj = m — rui — {h — 1). 


(3.3) 


Proof. By using the fact that 

deg adj_|_ = s — ruj, V s G {0,..., 2mj} 


and by comparing the principal degrees of both sides of (3.2). 

Lemma 3.3. If i,j G {1,..., n} satisfy rui + roj = h — 1 then c^J^' = 0. 

Proof. It suffices to consider the case i <n — 1. By Lemma 3.2 we have 


□ 





(3.4) 


0<mi+mj — {h—l)<2mj 

We will use the following orthogonality conditions [27] 


( 7*1 I ad^[^*2 7*=) = 0, if *17^*2- 


(3.5) 


It follows that 



0 . 


(3.6) 


On another hand, for ruj + ruj — (h — 1) > 2 we have 



(3.7) 
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For rrii + mj — {h — 1) = 1 we have 


{p'^ I ad/+ 7 -^) = (/+ I 7 ^) = 0 . (3.8) 

The last equality is due to i ^ n and so j cannot be 1, and /+ = const ■ 7 ^. Hence = 0. 

The lemma is proved. □ 


Proof of Theorem 1.6. Let M be any solution of (1.4). It can be uniquely decomposed in the 
form 

n 2mi 

M = ^5^/L,^adr^y (3.9) 

i=l m=0 

for some functions Kim = Kim{X). Substituting this expression into equation (1.4) we have 


n 2mi 


n 2mi 


Y, Y, K[„ adr,7‘+ E E +a 


2=1 m=0 


2=1 m=l 


n 2mi 


7 ' 


■, E E 


i=l m=0 


= 0 . 


(3.10) 


Here ' means the derivative w.r.t. A. Substituting (3.2) into (3.10) we hnd 

n 2mj n 2mj n 2mi 

E E(F.+^ E E E E ^im C-ijs 0 

j=l 5=0 j=l 5=0 i=l m=0 

where Kj _i := 0. Now using Lemma (3.2) we find 


(3,11) 


= -A". - A 


E 




''IJS 


K 


i,s+mi—mj-\-(h—l) 


mi>s—mj-\-(h—l) 


(3.12) 


where 0 < s < 2mj. Due to the inequality mt — mj + h — 1 > > 0 the system (3.12) is of 

triangular form. So we can solve it for the vector (.fT;>-i)J=i as long as the vectors (/T,s')j=i i — 
s + 1 are known. Note that, by dehnition Sj = Kj^ 2 m - Thus Kim, = 0,..., 2mj — 1 are 
determined as linear combinations of S'!,..., S'„ and their derivatives. Further noticing that 

s + mi — mj + {h — 1) > s + mi, \/i,j (3.13) 

we obtain the more subtle expressions (1.26). Indeed, it is easy to see that 

Kim = (-1)™ mi <m < 2mi - 1. (3.14) 

This observation together with (3.13), (3.12) yields (1.26). Note that the coefficients and 
k}muv ia (1-26) are uniquely specihed by the lowest weight structure constants of g. 

Taking s = 0 in Eqs. (3.12) we obtain a system of linear ODEs (1.27) of the form 

n 2mi n m-i 

d“’ = EE*^l'Sl’ + ^EE'=l.®l’' i = i.,,,.» (3,15) 

12=1 11=0 12=1 i ;=0 
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where the constant coefficients and are determined by the lowest weight structure con¬ 
stants of g. It is now clear that there is a one-to-one correspondence between solutions of the 
topological ODE (1.4) and solutions of the system (3.15). Thus these two systems of ODEs are 
equivalent. 

It remains to show that kl^Q = 0, G {1,... ,n}. The (s = 0)-equations in (3.12) read as 
follows 


K'io + A E 

mi-\-mj >h—l 


CjjO 


0 . 


Taking the A-derivative in the {m = 0, i = j)-equations of (1.26) we obtain 


(3.16) 


n 2mj-l n ruj-l 

— c(2™j + l) I ^ lO civ+1) Ul oiv) 

u=l i;=0 u=l i;=0 

n mj-l 

u=l v=0 


From equations (3.17), (3.16), (3.15) we know that Kjq does not contribute to kj^Q. The only 
possible terms that can contribute to kl^^^Q are 



mi—mj+(h—l) 

C-ijQ 


mi+rrij >h—l 


(3.18) 


Noticing again the inequalities rui — rrij + {h — 1) > rrii as well as (3.14), we find that (3.18) do 
not contribute to k\^Q unless 


rrii — rrij + {h — 1) = 2mj. (3.19) 

However, by Lemma 3.3 we know that = 0 if = h — 1. Hence kj^^ = 0. 

The theorem is proved. □ 

Equations (1.27) can be written in the following form (which will be used later) 

2 / 1-1 h 

Y, pI = 0’ ^ = (^1’ • • - Snf (3.20) 

li = 0 41=1 

where and p\ are constant n x n matrices, {p^)iu = k^uvi {pl)iu = k}^^. These matrices satisfy 
(P 2 mi+i) «, 2 mi+i = “1 and the following vanishing conditions 

(p°),„ = 0 if v>2mi + l, (3.21) 

{pl)iu = 0 if V >mi. 

Before proving Prop. 1.10 and Thm. 1.13 we do some preparations. Recall that the essential 
series Sa-,i of g are defined as the S'-coefficients of the generalized Airy resolvents M^. 
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Lemma 3.4. For any a G {1,..., n}, the essential series Sa-i, i = I,... ,n have the form 




C- G 


c 


m=0 


where Sai is a positive integer satisfying 

Sai = rua mod/i, Sai = rrii mod (h + l). 
In particular, the essential series Sa-a is non-zero and has the form 

OO 

5 ; c cS ^ 0, c e C, 

m=0 


(3.22) 


(3.23) 


(3.24) 


Proof. Fix an integer a G n}. By (1.17) we know that Sa-,i G A • C((A“^)), i = 

1,... ,n, where C((A“^)) denotes the set of formal Laurent series at A = 00 . If Sa-,i is zero then it 
has the form (3.22) with = 0, m > 0. Let Sa-i be non-zero. Recall the expansion (1.17) 


Ma{X) 


mg 

A 


Ma^k{X) 

k=l 


Ma,k{X) G L(g), degMa,fc(A) = m, 


{h + l)k, a = 1,..., n. 


Let Sa-,i,k{X) denote the S'-coefficients of Ma,k{X), we have 

OO 

Sa■,^ = A-"^ Sa-,,k{X). (3.25) 

k=0 

Noticing that deg ad^^'y* = m* we obtain 

deg Sa-,i,k{X) = ma-mi- {h + l)k, k>0. (3.26) 


Since deg X = h we conclude that 


ma — mi — {h -\- l)k = 0 (mod h). 


(3.27) 


Hence Sa-,i,k = 0 unless there exists an integer u such that k = ma — mi + uh, and in this case 

Sa-,i,k = const • ^3 28 ) 

which yields (3.23). Here, note that /c > 0 implies n > 0, hence Sai > 0. 

To prove (3.24) it suffices to show that Cg 7 ^ 0. It has been proved by Kostant that the 
elements a = 1,... ,n have the following decompositions^ (cf. p. 1014 in [27]) 

Am, = A/3nad^+“ 7 “, deg/3 =-(h - ma) (3.29) 

where v are some non-zero constant and /3 G g. Here, we have also applied the facts 

Ker ad/+ = Span{ad^+^ 7 ^,..., ad^™" 7 ”}, deg ad^™* 7 * = m*. (3.30) 

The lemma is proved. □ 

^We would like to thank Yassir Dinar for bringing this useful result of [27] to our attention. 
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Let us now prove Prop. 1.10. 

Proof. The statement 1) follows from (3.22). The statement 2) follows from (1.26) because Ma 
are linearly independent. The statement 3) is contained in Lemma 3.4. Let M be a regular 
solution to (1.4) such that the S'a-coefficient of M coincides with Sa-a- If the exponents of g are 
pairwise distinct then M must be equal to Ma since Sa-a is non-zero. Slight modihcations can 
be done for with n even as described in the content of the proposition. The proposition is 
proved. □ 

For any a G {1,..., n}, Thm. 1.6 reduces Ma consisting of n{h+ 1) scalar series to n scalar series 
J = 1, • • • Prop. 1.10 further reduces {Sa-j}f=i to one scalar series Sa-a- 

Finally let us consider the Laplace-Borel transform and prove Thm. 1.13. 

Proof. Equations (1.38), (1.39) follow directly from equations (1.26), (1.27). So i) is proved. 

To show ii), we will construct a basis of formal solutions to (1.39). Indeed, we consider the 
Laplace-Borel transform of the essential series Sa-i{X) 



(3.31) 


for suitable contours Ca- By Lemma 3.4 we have 



(3.32) 


where the coefficients are dehned by eq. (3.22). For i a, equations (3.23) imply that the 
integers Sai satisfy 


(3.33) 


Sai > rrta -|- h and Sai > rrii + h + 1 


if all exponents are pairwise distinct; for the case g is of the Dn= 2 k type, the above estimates are 
also true after possibly a linear change 


0n/2 0,1 0n/2 + ^2 <Pn/2+ly 

0n/2+l 0n/2 + ^2 0n/2+l 


(3.34) 

(3.35) 


of two solutions 0 „/ 2 , (pn/ 2 +i- III any case, for i = a we have 



(3.36) 


Hence we have obtained a set of independent formal solutions of of (1.39) 
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and we collect them to a matrix solution $ by defining ^ia = then we have 

0 \ 


/ X h ^ 0 

m 2 

^ = {I + 0{x)) ° " 

Vo 0 

Note that Equations (1.39) can be written as 

2/1-1 h 


0 


X h 


■-1 


x —)■ 0 . 


(3.37) 


rf </> + pI (j)' + vx'^ V) =0, (j)= (01,..., (pn? 


v=0 


v=l 


where and pi are n x n matrices defined in (3.20). So we have 

B{x) <l)'(a/) = C{x) <h(x), 


(3.38) 


(3.39) 


where B, C are n x n matrices whose entries are polynomials in x with coefficients determined 
uniquely by the lowest weight structure constants. More explicitly, the expressions of B, C are 
given by 

h 2h—l h 

c{x) = '^{-xY pl + ^{-iy vx'^-^pI. 

v=l f;=0 v=l 

We are going to show that the matrix B{x) is of the following form 


/ 0 ■ • ■ 0 bi ■ X \ 

: • •' b2 ■ X * 

0 .•■ .•■ ; 

\ bn ■ X * • • • * / 


(3.40) 


where the asterisks denote polynomials in x possessing at least a double zero at a: = 0. Recall 
that pI are defined as the matrix of coefficients of A x j = 1,..., n in the following equations 


K'jQ + A 


E 


"ijO 


^_2^7nj+7njH-/i—_ Q 


mi+mj'>h—l 


A 


(3.41) 


The summand with m* + = h — 1 in (3.41) does not contribute to the matrix B{x), since 

we have already proved in Thm. 1.6 that pi = 0. Let us consider the summands in (3.41) with 
rrii + rrij > h. Consider first the case g Y with n even where all exponents are distinct. In 
this case the {rrii + rrij = /i)-summand of the second term of l.h.s. of (3.41) contributes to the 
anti-diagonal entries of B(x). For mi + rrij = h + 1, h+ 2,... ,2h —2 the corresponding summands 
contribute to the entries of B (x) marked with an asterisk. The first term K'-q also contributes to 
B{x) through the last term of r.h.s. of (3.17), namely, 

n mj-l 

^ E E 

u=l 4;=0 
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(3.42) 


To understand this term we look at the (s = l)-equations of (3.12): 


^.0 = 



A 


E 

+mj 


nii-rrij+h, ^^rui+mi-h cimi+rrij-h) 
%1 


The contribution is of the form (3.40). We proceed in a similar way with analyzing Kji, Kj 2 , ..., 
in order to arrive at the matrix B (x) in a hnite number of steps. In the case of with n even, 
the above arguments remain valid with suitable choices of ^n/2^ ^n/2+l_ 

We have det 5 = ± 6i ■ ■ • x^. Since the dimension of the space of solutions of (1.4) is finite 

and since we have already obtained n linearly independent solutions of (3.39), we conclude that 
each bi is not zero. Using (3.21) and noting that the {rrii + mj = k, k > h + l)-off diagonal entries 
of B{x) (which are the terms marked with asterisk) have at least a {k — {h — l))-tuple zero, we 
conclude that the equations 

(j)'= B-^ C (j), 0 = (0i,...,0O^ (3.43) 


have the precise normal form (1.40), where the expression U_i is due to (3.37). 

The part ii) of Theorem 1.13 is proved. 

Finally, iii) follows from the normal form (1.40) automatically. The theorem is proved. □ 


The series (pa-i in the above proof are called the dual essential series of g. 

The hnal remark is about uniqueness of essential series. We £x a principal nilpotent element 
/+. Then, for a given choice of the basis { 7 *}(Li and of the series are 

determined uniquely. With a particular choice of only the element for any j G 

{ 1 ,... ,n} and hence each of the series is unique up to a constant factor except for 

the Dn with even n case. Because of this for presenting the expression of an essential or a 
dual essential series, we allow a free constant factor. For the exceptional case Dn with n even, 
7 ^/^, and hence {Si-n/ 2 -i *S'j;n/ 2 +i}r=i are unique up to invertible linear combinations with 

constant coefficients. 


4 Examples 

In this section we study examples of low ranks. We take a faithful matrix realization tt of g, and 
derive the normal form (1.40) and dual essential series. Denote by x the unique constant such 
that 


(a| 6 ) = xTr ( 7 r(a) 7 r( 6 )), Va, 6 Gg. (4.1) 

Below, we will often write 7 r(a), a G g just as a, for simplicity. 

4.1 Simply-laced cases 

Example 4.1 (The cases Ai,A 2 ). These two cases have already been studied in detail in Intro¬ 
duction. See (1.5)-(1.8), Example 1.9, Examples 1.15,1.16. Here we will make a short summary. 
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The essential series for Ai has the following explicit expression 


OO 


5=5: 

9=0 


(6^-1)!! 1 
969-gl 


which was derived in [4] by an alternative method. The dual fundamental series has the expression 


= X 2 e 12. 


The essential series for A 2 are given below in eguations (5.15)-(5.18). The dual fundamental 
series are given explicitly via confluent hypergeometric limit functions 



They satisfy the following dominant ODEs, respectively: 


27x2 - 81x0;.i + (x® - 84) 0i;i = 0, 
27 x 2 0".2 - 27 x 0 ' 2,2 + ( X ® - 21 ) 02;2 = 0 . 


Recall that the confluent hypergeometric limit function is defined by the convergent series 


oFi(;a;z) = ^ 
k=0 


1 

{a)k k\ 


where 


(a)fc = a(a + 1 )... (a + /c — 1) 

is the Pochhammer symbol. It is also closely related to the Bessel functions 

(f)' 


(4.2) 


Ju\^) - 


Ui' + i) 


O-Cl 


+ ) • 


Example 4.2 (The A^ case). In this case, g = sU{C), h = 4. The normal form of the dual 
topological ODE reads as follows 


-i 0 0 

0 -i 0 
0 0 -i 


+ 


X 


_A 

16 


4x2 



The three dominant ODEs read as follows 

64 (x^ - 15) x2 - 4 (125 - 3x®) x® - (x^® - 9x^° + 1134x® - 1260) fi-i = 0, 
4 X 02.2 + (x^ + 2) 02;2 = 0, 

64x2 0".3 - (64x - 12 X®) 03.3 - (x^° - 18x'^ + 36) = 0. 
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The dual fundamental series are given explicitly by 


01;1 = X 4 e 160 

, -1 
02;2 = X 2 e 20 , 




3.3 = X 4 e 160 . oFi 


4’ 4096 y 

3 \ 

4’ 40^7 


^ V 4’40%yJ ’ 


Example 4.3 (The A 4 case). In this case h = 5 and 0 = sl^lC). The normal form of the dual 
topological ODE of A^-type reads as follows 



0 

0 

0 \ 


( 0 

66 6 

25 ^ 

-^-fx 

\ 

0 

0 

3 

5 

0 

0 

2 

5 

0 

0 

X 

11 

210 

0 

0 

7 ™4 

300 

12 ™6 

5 

0 

0 

V 0 

0 

0 

-u 


\ 35 

0 

0 

0 / 


Let 9 = xdx'i then the dominant ODE for 04-4 reads 


(4.3) 


[3125 + 155) 9^ - 12500 {7x^‘^ + 620) 9^ + 625 {x^^ + 1277x^2 + 54870) 9^ 

-1250 + 321x^2 ^ 31124) Q + ^36 _ i495a;24 + 510995x^2 - 9215525)] 04;4 = 0. (4.4) 


Noting that x = 0 is a regular singularity of this ODE and that the indicial eguation at x = t] 
reads 



(4.5) 


one then recognises that the dual fundamental series is the (unigue up to a constant factor) 
solution corresponding to the root fc = — 


04;4 = X 



161 42 

_-I- 

210.35 


26 605 753 
223 .312.52 



In a similar way one can treat the other three dominant ODEs. 

Eor any solution 0 = (0i,... ,04)"^ of the normal form (4-3), 04 satisfies (4-4); Ihe series 

0a;4, a = 1,... ,4 are all non-zero. Hence, similarly to the proof of Prop. 1.10, for A 4 , one can use 
one dominant ODE to determine all solutions of (4-3). In particular, the interesting numbers 
Sa; 4 , a = 1, • • •, 4 (see eg. (3.22) above) can be read off immediately from (4-5): 


Si;4 — 16, S2;4 — 52, S3;4 — 28, S4;4 — 4. 


(4.6) 


We expect that this phenomenon always occurs when h is a prime number. 

Example 4.4 (The case Df). We have n = 4, h = 6 and mi = 1, m 2 = 3, m 3 = 3', m 4 = 5 and 

g = {B eMat{ 8 X)\B + SpB^pS = 0}, (4.7) 
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where rjij = Si+j^g, I < i, j < 8 and S = diag (1, —1,1, —1, —1,1, —1,1). In this case, X = ^ 


/0 1000000 \ 
0 0 1 0 0 0 0 0 
0 0 0 1 i 0 0 0 


A = 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 0 1 

I 0 0 0 0 0 0 

A 


0 0 
0 0 
1 0 
1 


yofoooooo/ 


The matrices Ai, A 3 , A 3 /, A 5 can he chosen as follows 


Ai = A, A5 = 2A^ A3 = A3, Ay = 2 


/0 0 0 0 1000 \ 

0 0 0 0 0 1 0 0 

0 0 0 0 0 0 1 0 

0 0 0 0 0 0 0 1 

AOOOOOOO 
OAOOOOOO 
OOAOOOOO 
\000A0000y 

The normal form of the dual topological ODEs of the D^-type reads as follows 



( 

6 

0 

0 

0 \ 


/ 26 6 
27 

0 

0 

91 ™3 1 1 ™10 

g X “T g a. 

\ 


0 

_ 1 

0 

0 

0 

0 

x^ 

0 

0 


0 ' = 

0 

2 

0 

_ 1 

2 

0 

-h 

X 

0 

0 

x^ 

0 



1 0 

0 

0 

-1 


V 

0 

0 

0 

/ 

The dual fundamental 

series 

are 

given by 







0 i;i = a: 6 exp 


13 

~Tm 


04;4 = X 6 



3 1 ’ 6 ’ 36 J ^ 63 6 ’ 36 


5 

6 ’^ 


They are solutions to the following dominant ODEs, respectively: 

108 (3a:^ + 182) - 4 (78 + 5461 + 9828) a: 

- (12 x^^ + 260 x^^ + 107029 a;^ + 62790) (fi-i = 0, 


^ 2;2 6 
+ X < 




3;3, 


V^ 2;2 - -r ^ V^ 2 ; 2 , V^3;3 “ — 2 ^ 

108 x^ 04.4 - (104 + 108 x) ^' 4.4 - (4 x^^ + 260 x^ + 39) 04;4 = 0. 


(4.8) 


(4.9) 


(4.10) 
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4.2 Non-simply-laced cases 

Example 4.5 (The case B 2 ). We use the realization given in [10]: 0 = {B G Mat(5,C) \B + 
SrjBrjS = 0}, x = | where rjij = 5i+j,6 (1 < j < 5), S = diag (1, —1,1, —1,1). We have 


A = 


/ 0 1 0 0 0 \ 

0 0 10 0 

0 0 0 1 0 

f 0 0 0 1 

V 0 A 0 0 0 / 


Recall that in this case h = 4, mi = 1, m2 = 3, and we can choose Ami follows 

Ai = A, A3 = 2AA 


The normal form of the dual topological ODE is given by 


-i 0 
0 -J 


X 


3 ^4 -7» J_ 2. 

4 ^ 16 ' 4 

0 


Solving this ODE system we obtain the dual fundamental series 


I —— - 

0i;i = a: 4 e 


Xi 


1 Scc'^ ^ 

02;2 = a;-4 ^ 


1 

’4’^ 
3 

’4’ 28 


5 X 


10 


+ 20 V’ 4’ 2" 


(4.11) 


(4.12) 


They are solutions to the following dominant ODEs 


(64x^ + 240) x^ 01.1 - (48x^ + 500) x® 0'i.i - (l6x^® + 36x^° + 1134x® + 315) 0i;i = 0, 

16 X^ 02-2 — 4 (3 X® + 4) X 02-2 — (4 X^® + 18 X® + O) 02;2 = 0. 

Example 4.6 (The B^, case). VTe have h = 6, g = {B G Mat(7, C) | B + SrjBrjS = 0}, x = 
where rjij = 1 < A j < 7 and S = diag (1, —1,1, —1,1, —1,1). The cyclic element is 


A 


/O 1 00000\ 
0 0 1 0 0 0 0 
0 0 0 1 0 0 0 
0 0 0 0 1 0 0 
0 0 0 0 0 1 0 
I 0 0 0 0 0 1 
\o fooooo/ 


ITe choose Ami follows Ai = A, A5 = 2 A®, A3 = \/2A^. The normal form of the dual 
topological ODEs read 


0 -i 
0 0 


0 

0 


M ™6 n 91 3 I 1 ™10 

rfiy a.' Vy a.' I r? 


X 


27 


2 2 
9 


X 

0 


0 

0 
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which is obviously equivalent to the normal form of restricted to 03 = 0. Hence 


0i;i = a; 6 exp 
02;2 = X 2 exp 
03;3 = X~^ exp 


13 
189 
x'^ 

13 

189 


0 1 1 > Q > 36 


13 


7x2: 0^1 ; 7; ^ 


63 


7 X 


14 


6’ 36 




5 x^^ 
’6’^ 


Example 4.7 (The G2 case). VTe have h = 6, mi = 1, m 2 = 5. Using the matrix realization 
given in [2], we obtain the normal form of the dual topological ODE of G 2 -type 


-- 0 

6 

0 -I 


-h 

X 


13 ™6 
54 

1 ^2 
360 


1820 ™3 I 5 ™10 

9 X ^ qX 
0 


(4.13) 


The dual fundamental series have the following explicit expressions 


_^ ±ox 

01;1 = X 6 e 756 


. 1. X^" 

0 1 > ’ 6 ’ 24 • 36 


, F (■ 

+ o-ci , 


252 


X 


14 


6’ 24 • 36 


02;2 = X 6 e 


13x* 

756 




X 


14 


6’ 24 • 36 


ITe remark that the normal form of G 2 -type coincides with the normal form of restricted to 
02 = 0 after a rescaling of K by A 2“3 A. 


Let us summarise in Table 2 the numbers ma^Va for the above examples. 


0 

h 

ma 

ra 

All 

2 

1 

1 

^2 

3 

1,2 

2,2 

^3 

4 

1,2,3 

2,1,2 

A 4 

5 

1,2, 3,4 

4,4,4,4 


6 

1,3, 3, 5 

2,1,1, 2 

B 2 

4 

1,3 

2,2 

B, 

6 

1,3,5 

2,1,2 

G 2 

6 

1,5 

2,2 


Table 2: Simple Lie algebras and intrinsic numbers. 


5 Concluding remarks 

5.1 On the analytic properties of solutions to the topological ODEs 

We have proved the dimension of the space of solutions regular at A = cxo of the topological 
equation of g-type 

[L,M]=0, L = dx + A (5.1) 
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is equal to the rank of g. We have also proved Thru. 1.6 on reduction of this equation through 
the lowest weight gauge. 

One can also consider the vector space of all solutions to the topological ODE; denote it by 
©(g). As A = 0 is a regular point for the system (1.4), every solution M(A) G ©(g) is uniquely 
determined by the initial data M(0). Thus dime ©(g) = dime g = n{h + 1). 

Denote the unique analytic solution of (1.4) determined by the following initial data 

M“^(0) = adj^ 7“, a = 1,..., n, A; = 0,..., 2ma. (5.2) 

Clearly M“^(A) form a basis of ©(g). With a particular choice of basis {Mi,..., M„}^ of ©^®(g), 
there exist unique partial connection numbers Cabk such that 

n 2mf, 

= (5.3) 

b=l k=0 

where Cabk are constants. We will study these connection numbers for the topological ODEs as 
well as the monodromy data for the dual topological ODEs in a subsequent publication. 

5.2 Application of topological ODEs to computation of intersection 
numbers on the moduli spaces of r-spin structures 

First, it should be noted that all main results of the previous sections remain valid after a change 
of normalization of the topological equation. It will be convenient to introduce a normalization 
factor K in the following way 


[L,M]=0, L = dx + KA. (5.4) 

This parameter will be useful for applications of solutions to (5.4) to computation of certain 
topological invariants of Deligne-Mumford moduli spaces. In particular, here we explain two 
applications of the topological ODEs (5.4). The proofs follow the scheme used in [4] for computing 
the coefficients of Taylor expansion of the logarithm of the Witten-Kontsevich tau-function of 
the KdV hierarchy. Details of the proofs can be found in an upcoming publication [5]. 

Let M.g,N denote the Deligne-Mumford moduli space of stable curves C of genus g with 
N marked points xi, ..., xn, jCk the /c*^-tautological line bundle over i^k = Ci i^k) ^ 

Let r > 2 be an integer. Fix a set of integers {ak}k=i C {0,..., r—1} such that 2,g—2—'Y^^=i a-k 
is divisible by r. Then the degree of the line bundle Kc — divisible by r. Here Kq 

is the canonical class of the curve C. Hence there exists a line bundle T such that 

N 

c^Kc-J^aiX,. (5.5) 

i=l 

^Such a choice depends on a choice of an appropriate sector in the complex plane near A = 00 . 
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Such a line bundle is not unique; there are choices of T. A choice of such a line bundle T is 
called an r-spin structure on the curve {C, Xi,... ,xn) G Aig^w- The space of all r-spin structures 

/ T 

admits a natural compactihcation Aig-^ai,...,aN with a forgetful map 

V : Adiai,..,ajv Aig^N- 

- 1 It 

There is a natural complex vector bundle V —?■ -Alg;ai,...,ajv hber over the generic 

point {C,xi,... ,xn,T) G -Alg;ai,..,,ajv coincides with H^{C,T). Here the genericity assumption 
means that H^{C,T) = 0. Under this assumption, using the Riemann-Roch theorem one con¬ 
cludes that rank V = s + g — 1 where 

czi . "T cztv — (25^ — 2) 
r 

The construction of the vector bundle V over non-generic points in AAg.^^ is more subtle; for 
details see [16] and references therein. The Witten class cw{cbi, ■ ■ ■ ,a,N) £ is 

dehned via the push-forward of the Euler class of the dual vector bundle 

cw{ai ,..., Oat) = (euler (V^)). (5.6) 

'f9 

See more details on properties of cw{ai ,..., ajy) in [38, 37, 16, 22, 32]. 

Witten’s r-spin intersection numbers are nonnegative rational numbers dehned by 


{riiM ■ ■ ■ '^iN,ki,Y ■= [_ A-■■ A Acw{ii-l,---AN - !)• (5.7) 

•^g,N 

Here ii G {1, ... ,r — 1}, ki > 0, i = 1,..., N. Note that we have assumed that the genus in l.h.s. 
of (5.7) is reconstructed by the dimension counting: 

——- + -h ——- + -— -{g - 1) + ki-\ - \- kjsf = 3g - 3 + N. (5.8) 

f TV 

See [6, 7, 29, 28, 22, 39, 3, 8] for some explicit calculations of Witten’s r-spin intersection numbers. 

According to Witten’s r-spin conjecture [38] the intersection numbers (5.7) are coefficients 
of Taylor expansion of logarithm of tan-function of a particular solution to the Drinfeld-Sokolov 
integrable hierarchy of the An type, n = r — 1. This conjecture was proved by Faber-Shadrin- 
Zvonkine [16]. 

We will now use the constructed above generalized Airy resolvents for a simple algorithm for 
computing the r-spin intersection numbers in all genera. 

For a given A^ > 1 and a given collection of indices R, ..., satisfying 1 < < r — 1 dehne 

the following generating function of the Wpoint Witten’s r-spin intersection numbers by 


..., Ajv) = (-K 


N 


E ( 


-1) 


+ + 


N (ii) 




(T„.,...r.,,.,)’'-*n5.9) 


Here Ai, ..., Aat are indeterminates, k = ^■ We use the standard notation (4.2) for the 

Pochhammer symbol. 
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Theorem 5.1. Let n = r — 1, g = s/„+i(C), A = Y^=i Ei,i+i + XEn+i,\- Let Mi = Mi{X) be the 
basis of generalized Airy resolvents of Q-type, uniquely determined by the topological ODE 

M' = k[M,A], (5.10) 


normalized by 


Mi = —A Aj + lower degree terms w.r.t. deg, Aj := A*. 


(5.11) 


Denote rjij = 6i+j^n+i- Then the generating functions (5.9) of N-point correlators ofr-spin inter¬ 
section numbers have the following expressions 


dF: 


r—spin 


7' — 1 

■ ('^) ^ (-'^2) 1,71+1 ('^) X ^ *^7,715 


dA 

T^r-spinf^ y \ _ _+ Tr {Mj^^ (AgJ ■ ■ ■ Mj^^ {Xsn)) 

1) • • • ) N) J^r 2^ \ \ 

sgSn 


(5.12) 


-dN2 Viii2 


Ai '* A2 (A Ai + *2 A2) 

(Al - A2)2 ^ 


N>2. (5.13) 


Note that up to a constant factor the entry (Mj)i „+i is nothing but the essential series Si-n- 
Then, from (5.12) we expect that all one-point r-spin intersection numbers can be read off from 
coefficients of solutions to the dominant ODE expanded near the regular singularity a; = 0. 
We also remark that alternative closed expressions for one-point r-spin intersection numbers have 
been obtained in [28] by using the Gelfand-Dickey pseudo-differential operators. 

Example 5.2 (r = 2). Witten’s 2-spin invariants coincide with intersection numbers of fj-classes 
over Aig^N [36, 26, 16], So Thm. 5.1 in the choice r = 2 recovers the result of [4, fO]-' for N >2, 


Y, Y1 (2pi + 1)!!.. ■ (2p„ + 1)!! 

3=0 pi,...,pjv>0 


2 pi +3 2 pjY +3 

.PI . . . ol,PN \ 2 . . . y 2 


N 


-- y. 

N ^ 


J Mg^N 

Tr(M(A,J---M(A,^)) 


Al ^A2 ^ (Al -|- A2) 


P.5. nAi(A^ 


r,. ■ Kj+x, 


-6 


N2- 


where M = 


1 


/ (63-5)!! \-33+2 

Z^3=l (96 k2)®"^-(9-1)! 


(Al - A2)2 

9 (6g-l)!! A-3q \ 

2 -^g={) ( 96 K[^) 3 -g\ 


(5.14) 


k = (-2) 


-1 


nS^^oo 63+1 (63-1)!! \-3o+l _ (63-5)!! \-3o+2 

\ W / g =0 63—1 (96 k^)9'3! 2_/3=1 (96 k2)9“1-(3—1)! / 

The formula (5.14) is closely related to a new formula recently proved by Jian Zhou [ 4 I] (See 
Thm. 6.1 therein). For N = 1, it follows easily from (5.12) the well-known formula 


('^33-2)3 — 


249 ■ g\ 


for g>l. 
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i 

( 2Kii + K 22 

- 2^1 - 37723 

6 .S 2 

(r = 3). Define M = 

277io + 7721 

—2/722 

37723 - 2 .S 1 

1 

V 7720 

1 

0 

CM 

—277ii + K 2 

77io = ^ + 3A^2, 7^11 = --, 

K 



o(4) 

K2o = ^-2XSi, 

cw 

-^21 —- 

Off 

J^22 — — 

0 / 

7^23 —- • 


The generalized Airy resolvents Mi (A) and M 2 (A) are obtained by replacing the functions Si and 
S 2 in the above expressions by 


Sl;l = 


Sl;2 = 


and 


S2;l = 


1 (-i)^-37^r(8(7 + i) 24..+1 

2233^2 ^ 249 ■ g \ ■ ( 54 fi: 2 ) 3 ff-i r(^ + 1 ) 

1 ^ (-l)^-3^^r(8f7 + f) 242^ 

2332^2 ^ 249 . g\ ■ (54 T{g + |) 


^ ^ (-l)9.3"9r(8^+f) ^_24.+ir 

23 ^ 249.^!.(54«:2f"+2p^^^|) ^ ’ 


9=0 ^ - vy ' 3^ 

00 / 1 o7n -n/o _ , 2' 


^2;2 - 

9=0 


-i)9-3^9r(8^ + 


6 j^^249 ■g\-{5AK^)^9 .T{g + l) 
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respectively. Setting k = ^ we easily obtain from Theorem 5.1 that 

1296 


( 9 - 1 , 89 - 7)9 - 


('r2,89-2)9 - 


466569-(^-I)! - (I) 
1 


{9 > 1), 


466569 ■ (1)^ 

1 

N 


(9 > 1 ), 


^ \ __^ Tr (AsJ ■ ■ ■ (^ sn )) 

1) • • • ) jy- 2_^ T-rAf / \ \ \ 

seSN 


(5.15) 

(5.16) 


(5.17) 

(5.18) 


Ai A2 (*1 Ai + 12 A2) 

-0N2 Viii2 -7T-, N >2. 

(Ai — A2)^ 

The above formulae of one-point 3-spin intersection numbers agree with known formulae in [Sf, 
1, 29, 28]; e.g. the first several of these numbers are given by 

(ri,i)i = —, (r2,6)i - ^Yio4’ “ 746496 ’ ^^2,14)2 - 4337294080’ 




















5.3 Topological ODEs, Drinfeld—Sokolov hierarchies, and FJRW “quan¬ 
tum singularity theory” 


According to [17, 30], the partition fnnction of FJRW invariants for an ADE, singularity 
with the maximal diagonal symmetry group is a tan-function of its mirror Drinfeld-Sokolov 
hierarchy; the partition function of FJRW invariants for a {D 2 k, (J)) singularity is a particular 
tan-function of the iJ 2 fc-E)rinfeld-Sokolov hierarchy; the partition function of non-simply-laced 
analogue of the FJRW intersection numbers is dehned as a tan-function of the Drinfeld-Sokolov 
hierarchies of B C F G-tjpe [30]. 

Definition 5.4. The Drinfeld-Sokolov partition function Z is a tau-function of the Drinfeld- 
Sokolov hierarchy of Q-type [10] uniquely specified by the following string equation: 


n oo 




dZ 1 ^ i ,■ 




dtl' 


(5.19) 


Here, tl, are time variables of the Drinfeld-Sokolov hierarchy and rjij = 5i+j^n+i- 


It should be noted that the terminology “the Drinfeld-Sokolov hierarchy of g-type” that we 
use refers to the Drinfeld-Sokolov hierarchy (under the choice of a principal nilpotent element) 
associated to the non-twisted affine Lie algebra Since the construction of this integrable 
hierarchy does not depend on the central extension of the loop algebra, it is essentially associated 
with the simple Lie algebra. 

Dehne the following generating functions of iV-point correlators of Z by 


.,IN 




(0)(5.20) 


fcjv 


We will now express these generating functions in terms of the generalized Airy resolvents of 
g-type. To this end we will need to use the following multilinear forms on the Lie algebra 

R(ai,... ,aAr) = tr (adai o ... o ada^), Voi,..., oat e g. (5.21) 

Theorem 5.5. Let g be a simple Lie algebra of rank n. Let Mi = Mi{\), i = 1,... ,n be the 
generalized Airy resolvents of g-type, which are the unique solutions to 

M' = n[M,A], K=(yEh^ ^ (5.22) 


subjected to 


Mi = —A h Arm + lower degree terms w.r.t. deg. 


(5.23) 


Then the generating functions (5.20) for the N-point correlators of the Drinfeld-Sokolov partition 
function associated to g are given by the following expressions 


dA 


(A) 


-^B{E_e,Mi) - kX 6i^n 


(5.24) 
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• • • 5 


E 




2NK^ ^ rr^ fA - A ) 

^il '^12 

A^ '* A2 '* (mil -^1 + "^*2 -^2) 


-(5Ar2 flili2' 


(Ai-A2)2 


iV > 2. (5.25) 


i/ere, hy is the dual Coxeter number. 


Similarly as before, we expect that all one-point correlators of the Drinfeld-Sokolov partition 
function can be read off from coefficients of solntions to the n-th dominant ODE expanded at 
X = 0. The proof will be given in [5]. It is interesting to mention that for the ADE cases, 
the correction term appearing in (5.25) for N = 2 coincides with the propagators derived in 
[31] in the vertex algebra approach [1] to FJRW invariants; in these cases, the Drinfeld-Sokolov 
partition fnnction coincides with the total descendant potential [19, 15, 14] of the corresponding 
Frobenius manifolds [11, 12]. It would also be interesting to investigate relations between the 
dual fundamental series and the hypergeometric functions considered in [21]. 

In subsequent publications we will continue the study of topological ODEs as well as their 
applications to computation of invariants of the Deligne-Mumford moduli spaces considering 
generalized Drinfeld-Sokolov hierarchies [20]. 


Appendix A Simple Lie algebras and generalized Airy 

functions 

Let 0 be an arbitrary simple Lie algebra of rank n. We consider in this appendix a g-generalization 
of equation (1.10). 

Let E_g and /+ be dehned as in the Introduction. Let A \E_g. Throughout this 

appendix we will £x a matrix realization vr : g —)■ Mat (m, C). Again for any A G g we simply 
write 7r(A) as A. 

Definition A.l. The system of m linear ODEs 

^ + Ay = 0, ( 2 / 1 ,... , 2 /™)^ e C"*, A e C (A.l) 

dX 

is called the generalized Airy system of Q-type. 

The dimension of the space of solutions to (A.l) is m. Let F : C —)■ Mat (m, C) be a fundamental 
solution matrix of (A.l), i.e. 

W + A F = 0, det (F) = const ^ 0. (A.2) 

Proposition A.2. Let A be any constant matrix in vr(g). The matrix-valued function M := 
FAF“^ is a solution of the topological ODE, i.e. we have 

M' = [M,A]. (A. 3) 
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Proof. By straightforward calculations. 


□ 


Noting that A = 0 is a regular point of the generalized Airy system (A.l), we dehne a 
particular fundamental solution matrix Y of (A.l) by using the following initial data 

1^(0) = Im (A.4) 

where Im denote the m x m identity matrix. Then we have 

Proposition A. 3. Let {M‘^^}a=l^...,n,k= 0 ,..., 2 ma denote the basis of &{g) defined by (5.2). Then 

= Y 7r(ad^^7“) (A.5) 

Proof. By using Prop. A.2 and by using the standard uniqueness theorem for linear ODEs. □ 
Corollary A.4. Let {Ma}a=i,...,n be a chosen (analytic) basis of Then we have 


71 


n 2mi, 

b=l k=0 


(A.6) 


where Cabk o,re the corresponding partial connection numbers (c.f. (5.3)). 


Example A.5 (A„, n > 1). Take g = s/„+i(C). We have m = n + 1. Any solution y = {yi = 
y,y 2 ,..., ymf of (A.l) satisfies that 


(-l)"+i|/("+i) = A|/, 

= k = 2,...,m. 

Solutions to (A.7) can be represented by a Pearcey-type integral 

f.n+2 


y{X) = / exp 


X 


\x] dx 


(A.7) 

(A.8) 


(A.9) 


Ic \n + 2 
where C is a suitable contour on the complex x-plane. 

Example A.6 (iA„, n > 4). Take the matrix realization of q as in [10]. We have m = 2n. The 
cyclic element A takes the form 


/ 0 1 


A = 


\ 


A 

2 , 
VO I 


1 

0 / 
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For any solution y = {yi = y, y 2 , ■ ■ ■, ymY ^o (A.l) we have 

(A.10) 

yk = {-lf~^y^^~^\ fc = 2,...,n-l, 

(All) 

Vk = (-1)''“^ k = n + 2 ,...,2n-l, 

(A. 12) 

2 v'„ = »' +. = (-1)”-' -fiy, 

(A.13) 

Solutions to (A. 10) can be represented by a Pearcey-type integral 

9 (A) exp ( 

(A.14) 

where C is a suitable contour in the complex x-plane. 

Example A.7 {Bn, n >2). Take the matrix realization of q as in [10]. 

We have m = 2 ?7, + 1, 


and g = {B G Mat(m, C) | B + SrjBrjS = 0} where rjij = 5i+j,2n+2 (1 < hj < 2n + 1), S' = 
diag (1, —1,1, —1,..., 1, —1). We have 



/ 0 

1 0 

0 \ 

A = 

0 

0 ■■■ 

■■■ 0 


A 

2 

0 ■■■ 

■ ■■ 1 


V 0 

A 0 

2 ^ 

0/ 


Then any solution y = {yi = y, y 2 , ■ ■ ■, ymY' of (A.l) satisfies 


^(2-+i) = A2/' + -2/, 

Vk = (- 1 )^"^ k = 2,...,m 

X 


1 , 




y- 


(A.15) 

(A.16) 

(A.17) 


Solutions to (A.15) can be represented by the following Pearcey-type integral 

4n+2 


y{X) = / exp 

Jc 


X 


2n + 1 

where C is a suitable contour in the complex x-plane. 


\x ] dx 


(A.18) 


Example A.8 (Eg). We take the matrix realization the same as in [13], where m = 27. The 
matrices J+ and E_g can also be read off from [13], The corresponding generalized Airy system 
reduces to two linear ODEs for yi, y^. Denote by u{x), v{x) the Laplace-Borel transforms of 
yi, i/g respectively: 


yi{X) = / u{x)e ^"^dx, 2/6(A) = / v{x)e 

Jc Jc 


—X X 


dx. 


(A.19) 
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Then u, v satisfy the following system of ODEs 


u 


It follows that 



0 


2 

3 



5x^^ 

3 

26 „16 

9 

5x^^ 



-27 u"{x) + (270 + 189 x) u' + + 675 x^^ + 75)u = 0, 

-27x"(x) + (270 x^^ + 405 x) v' + (x^® - 405 x^^ + 300) x = 0. 


Solving (A.21) we have 


_1 5^13 

u = cix 3 e 13 0-^1 



25 5 3:13 

+ C2 X 3 e 13 0-^1 



where ci,C 2 are arbitrary constants. 


(A.20) 


(A.21) 

(A.22) 


(A.23) 


Applying the saddle point technique to the above integral representations one can derive 
asymptotic expansions for the generalized Airy functions. In principle, Proposition A.2 can be 
used for computing the asymptotic expansions of solutions to the topological ODE. However, the 
methods explained in the main part of the present paper seem to be more efficient for obtaining 
such expansions for regular solutions. 

Remark A.9. For the original construction [26] of the Witten-Kontsevich tau-function it was 
introduced a matrix analogue of the Airy integral. More recently in [18] the Kontsevich’s matrix 
Airy function was generalized for an arbitrary compact Lie group. For g = s/„+i(C) the gener¬ 
alized matrix Airy function appears in [9] in representation of the tau-function of the Drinfeld- 
Sokolov hierarchy associated with the r-spin intersection numbers on Aig^N, n = r —1 (see above). 
It is not clear whether the generalized matrix Airy functions of [18] for other simple Lie algebras 
can be used for computing the FJRW intersection numbers. 
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